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We show that doping a Majumdar-Ghosh chain with non-magnetic impurities does not produce 
almost free spins 1/2. The difference between this system and other spin liquids, such as unfrustrated 
spin ladders, is illustrated in the context of the general dimerized, frustrated spin chain. By detailed 
analysis of the excitation spectra of finite chains with two impurities, we investigate the evolution of 
the screening affecting impurity-induced free spins, and of the sign and magnitude of their effective 
interactions. We comment on a possible connection to impurity-doping experiments in CuGeOa. 



PACS numbers: 75.10.Jm, 75.30.Hx, 75.40.Mg, 75.50.Ee 
I. INTRODUCTION 



Doping of low-dimensional spin systems by static, non- 
magnetic impurities has been a topic of increasing inter- 
est in recent years, driven by the observation that signifi- 
cant information may be obtained concerning the nature 
of magnetic fluctuations. For antiferromagnets- the issue 
was first addressed in the Heisenberg modeljil where a 
small induced moment was found around each empty site, 
A comprehensive study of Hubbard antiferromagnet: 
suggested that the size of the effective moment may vary 
from this small value in the limit of large U/t to a full 
compensation of the missing spin at small U/t. Nonmag- 
netic impurities in high-temperature superconductors be- 
yond the antiferromagnetically ordered region at low hole 
doping were found in nuclear magnetic resonance (NMRA 
measurements to give localized moments of S = 1/2, El 
while further NMR studies have investigated the effects 
of impurities on magnetic correlations^ and shown that 
each vacancy causes a strong antifermmagnetic (AF) po- 
larization of its immediate vicinity.El A full theoretical 
discussion can be found in Ref. ^. 

Nonmagnetic impurities in quantum spin liquids give 
rise to various phenomena, among which early attention 
was focused on the formation at low temperaturexifJarge- 
spin clusters with random, effective interactionaJu on a 
subsystem of fcee spins isolated by doping. For ladders 
and dimerizedlj geometries these interactions are coher- 
ent and quasi-long-ranged, in accord with the weak anti- 
ferromagnetism observed in the low-temperature phases 
obtained pn Zn-doping of the quasi-ld ladder material 
SrCu203,til and of the dimerized chain CuGe03.t3 These 
ideas have been extended to discuss true magnetic or- 
der on the impurity spin subsystem in higher dimensions 
for coupled, dkacrized chainstJ and for the 1/5-depleted 
square lattice,Ej both_pf which are bipartite. More gen- 
erally, Martins et al&3 argued for local enhancement of 
AF correlations around impurity sites in a range of low- 
dimensional S = 1/2 systems, and later discussed the 
evidence for nonmagnetic impurities liberating effectively 



free spins ,0 a notion which now has broad acceptance for 
all spin liquids. 
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FIG. 1. The dimerized, frustrated spin chain. 

However, very few of these studies have considered sys- 
tems with magnetic frustration. In this paper, we show 
that in frustrated systems missing moments do not neces- 
sarily produce nearly-free spins 1/2, even when the spin 
gap is large. The analysis is based on the Heisenberg 
model for a chain with dimerized nearest-neighbor and 
frustrating next-nearest-neighbor interactions (Fig. 1), 
described by the Hamiltonian 

H = Jr J2 [(1 + 8)S*-S* + (1 - <y)S?.S? +1 



+j 2 s ?- s ?+v 

i,i;m— a,b 



(i) 



We perform zero-temperature Lanczos diagonalization of 
finite S — 1/2 chains containing two S = impurities, 
and from inspection of the low-energy excitation spec- 
trum deduce the nature of the impurity-induced ground 
state. All of the results to follow were obtained on pe- 
riodic chains of 20 sites, meaning with 18 spins and 2 
impurities, but we have verified with longer chains that 
the results on which we will base our conclusions are dom- 
inated by the shorter inter-impurity separation, and re- 
main essentially independent of the chain length. 



1 




FIG. 2. Schematic representation of nearest-neighbor 
dimer singlet coverings in the ladder (a) and in the frustrated 
chain (b,c) with two impurities, represented by the open cir- 
cles. 



II. SPIN LADDER 

We begin by illustrating the nature of our analysis for 
the well-known case of the spin ladder, which is the 6=1 
limit of Eq. (|]J). As shown in Fig. 2(a), each impurity 
is expected to leave an effectively free spin on the de- 
pleted rung. By "effectively free spin" we understand an 
5=1/2 degree of freedom whose interactions with any 
other moments in the system occur on an energy scale 
very much smaller than the superexchange interactions 
in Eq. (|l|). In a ladder the interactions between the free 
spins are exponentially weak because of the spin gap and 
consequent short correlation length, and have an effective 
AF or ferromagaetic (FM) nature according to the Mar- 
shall sign ruleJlZMa In Fig. 3(a) are shown the lowest four 
energy levels for ladders with all chain-to-rung coupling 
ratios from to 1, for a configuration with the first impu- 
rity at (l,b) (first rung, b-chain; see Fig. 1) and the sec- 
ond at (4,b). Because these are on opposite sublattices, 
the lowest level is a singlet, while the first excited state is 
a triplet, and there is always a wide separation of this low- 
est manifold from other levels. This situation is mirrored 
for the impurity configuration ((l,b)(4,a)) [Fig. 3(b)], but 
with singlet and triplet levels interchanged. We interpret 
these results as corresponding to two free spins isolated 
by the impurities, and communicating with effective AF 
(a) or FM (b) interactions given by the singlet-triplet 
separation in the lowest manifold. 

In Fig. 4 we show the singlet-triplet gap A$t as a func- 
tion of impurity separation for ladders with J2/2J1 = 
0.25 and 1. The alternation of the sign of A$t is clear, 
as is the decay in its magnitude with separation. For 
J2/2J1 = 0.25 this decay is very rapid, in keeping with 
the expectation of a short coherence length on an un- 
doped ladder with a large singlet-triplet gap. We note 
that the gap values with impurities on adjacent rungs 
correspond approximately to two sites with effective cou- 
pling J2 when they occupy the same chain, but be- 
come very small when they are on opposite chains as 
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FIG. 3. Lowest four energy levels Ei as a function of the 
chain to rung coupling ratio, J c /J r , for the spin ladder with 
impurities on the first and fourth rungs. In (a) the impurities 
occupy the same chain and in (b) they are located on oppo- 
site chains. The labels denote first and second singlet and 
triplet states. For the ladder defined by Eq. ([!]) with 5 = 1, 
Jc/Jr = J2/2J1. 



in this case the ladder is cut. For the isotropic ladder, 
J2/2J1 = 1, the decay is no longer clearly exponential, 
and implies that despite the spin gap in this system the 
impurity-induced states have a significant spatial extent 
which is approaching the system size. Our studies of the 
energy spectra provide a clear and complete confirma- 
tion of the results of Sigrist and FurusakH for the sign 
and magnitude of effective interactions between free spins 
in the ladder. 



III. FRUSTRATED CHAIN 

We turn now to the frustrated chain, and begin with 
the case of no dimerization (6 = 0). Much is known 
about this system, including the presence of a quantum 
critical point separating gapless and gapped phases at 
J2/J1 = 0.2412, and the existence of an exact wavefunc- 
tion at the Majumdar- Ghosh (MG) point J2/J1 = 1/2. 
The gapped phase arises due to spontaneous dimeriza- 
tion, and possesses a double degeneracy on translation. It 
is immediately apparent from Figs. 2(b,c) that the degen- 
erate dimer patterns can simply adapt to accommodate 
a single impurity, while two impurities may fall in one of 
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impurity rung 

FIG. 4. Singlet-triplet energy gap Ast for ladders with two 
impurities. The first impurity is located at (l,b), and the sec- 
ond on the rung indicated on the z-axis, on the same (dia- 
monds, triangles) or the opposite (squares, circles) chain. 

two ways. If their separation is even, the odd number of 
sites between them has no ideal covering (one dimer must 
always fall on a J 2 bond), and two nearly-degenerate sin- 
glet states are expected [represented by the solid and 
dashed dimers in Fig. 2(b)]. If the separation is odd, the 
even number of sites between the impurities can be ide- 
ally covered by the solid dimers in Fig. 2(c), whereas the 
second singlet state represented by the dashed dimers is 
driven to higher energy. 

Fig. 5(a) shows the lowest levels in the energy spectrum 
for the case of impurities located on sites (l,b) and (6,b), 
a situation analogous to Fig. 2(b). For J2/J1 > 1/2 we 
do indeed find a pair of low-lying singlets as expected 
from the above argument, and in fact this manifold lies 
well below the closest triplet states for all values of the 
frustration away from J 2 /Ji = 1/2. However, below the 
MG point there is an abrupt level crossing such that the 
lowest manifold consists of a singlet and a triplet. These 
become degenerate at J 2 = 0, which corresponds to the 
limit of the doped Heisenberg chain: in this case the 
chain is simply decoupled into two odd- length segments, 
each an effective S = 1/2 degree of freedom, while if the 
impurity separation is odd (below) the two even-length 
segments are 5 = 0. For a long chain randomly doped 
with a finite concentration of impurities.— one observes 
on average half a free spin per impurityJIS The results 
for frustration values of < J2/J1 < 0.45 imply that the 
chain segments lose the resonant character which permits 
two low-lying singlets, resulting in a pair of free spins 
with a single effective coupling which decreases linearly 
with J 2 . 

Fig. 5(b) shows the spectrum for the case of impurities 
with odd spacing, where one singlet state lies well below 
all of the other levels over the entire range of frustra- 
tion, also as expected from the above considerations. For 
appreciable frustration this is the favored (solid) singlet 
state of Fig. 2(c), while for weak J 2 it has more the char- 
acter of two separated, short chains of even length. In 
no circumstances do the impurities introduce free spins. 
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FIG. 5. Lowest four energy levels Ei as a function of frus- 
tration J2/J1 for the undimerized chain with two impurities 
whose separation is even (a) or odd (b). The labels denote 
first and second singlet and triplet states. 

Thus for a random distribution of impurities, strongly 
frustrated chains will not show free spin degrees of free- 
dom, but would exhibit a crossover to half a free spin per 
impurity as frustration is reduced. 

In the light of our opening statement concerning size 
effects in our calculations, it is clear that impurity doping 
destroys the quantum critical point at J2/J1 = 0.2412, 
which might have been expected to mark the boundary 
between these two types of behavior, and replaces it with 
a crossover. At very low impurity concentrations, this 
critical value should still mark the onset of the crossover 
as a function of J2/J1. By contrast, for the high doping 
concentrations represented by our results, the effect of J 2 
in localizing a spin of 1/2 at either end of a chain segment 
cannot compete with the finite-size gaps in the spectra of 
short segments, the regime of half a free spin per impurity 
is moved to larger J 2 ~ 0.5, and the crossover broadened. 

Before leaving the undimerized chain, we present 
briefly the effects of altering the impurity spacing in the 
MG chain. Fig. 6 displays two distinctive features in 
sharp contrast to the case of the spin ladder (Fig. 4). 
The first is that the alternation of sign in the ladder case 
is no longer present, being replaced by the alternation 
in character of the ground state (one favored singlet, or 
two similarly disadvantaged singlets) with odd or even 
impurity spacing. The second is that the energy level 
separations have no appreciable dependence on impurity 
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4 5 
impurity rung 

FIG. 6. Singlet-triplet energy gap Ast for MG chains 
(J2/J1 = 0.5) with two impurities. The first impurity is lo- 
cated at (l,b), and the second on the rung indicated on the 
x-axis, separated by an even (squares) or an odd (triangles) 
number of sites. Also shown (circles) is the energy Ass of the 
first excited singlet state, which is nearly degenerate with the 
ground state for even impurity separation. 



spacing beyond this alternation, demonstrating clearly 
the local dimer nature of the states involved. We note 
in passing that while nonmagnetic impurities in a spin 
ladder destroy its spin-liquid nature by creation of states 
within the spin gap associated with the free spins, the 
dimer reorientation in the MG chain effectively acts to 
preserve the spin-liquid character in the sense of retaining 
a significant singlet-triplet gap. 



IV. DIMERIZED, FRUSTRATED CHAIN 

We turn now to the general case of the dimerized, frus- 
trated chain, with which one would wish to explore the 
change in behavior from the MG regime to the ladder 
regime. We investigate this change by varying the dimcr- 
ization parameter 8 from to 1 at a fixed value of J 2 . Wc 
will illustrate our results with the case Ji — 0.5, which 
connects the MG chain to a ladder of chain-to-rung cou- 
pling 1/4, and have verified from studies elsewhere in 
(J2, 8) parameter space that the behavior shown is quite 
generic. Fig. 7 gives the low-energy spectra for a selec- 
tion of impurity configurations, presented in the form of 
the energy level separations Ast and Ass between the 
lowest singlet state, which gives the zero of energy, and 
the first triplet and second singlet, respectively. The sec- 
ond triplet level, which is not shown in Fig. 7, lies close 
to the second singlet only at 8 = (Fig. 5). 

In Figs. 7(a) and (b) the gaps are shown for impurities 
on the first and fourth rungs, and in Fig. 7(c) for the first 
and fifth. At small dimerization Ast is large, and for im- 
purities on the same chain the first excited singlet lies in- 
side this gap. However, this level is driven to much higher 
energies by even rather small dimerization, and for most 
of the parameter range the low-energy manifold contains 
only one singlet and one triplet. The sharp increase in 
Ass is readily understood from the fact that dimeriza- 
tion of the Ji bonds limits the flexibility of the system to 



to 

+ 



5) 



0.50 
0.40 
0.30 
0.20 
0.10 
0.00 
-0.10 



/ 
/ 

/ 


A ST ((1,b)(4,b)) 

A ST ((1,b)(4,a)) 

A ss ((1,b)(4,b)) . 






(a) 



0.0 0.2 0.4 0.6 0.8 1.0 

8 



to 
+ 



CM 



0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 
-0.03 



1 \ 1 

V \ _ 

V \ 


A ST ((1,b)(4,b)) 


\ \ _ 
^ \ 


A ST ((1,b)(4,a)) ^^^^ 


\ \ 

M 








" (b) 







0.0 



0.2 



0.4 



0.6 0.8 

S 



1.0 



to 
+ 



5a 



0.010 



0.005 



0.000 



-0.005 



a st «i,b)(5,b)) 
a st «i,b)(5,a)) 




FIG. 7. Singlet-triplet energy-level separation, Ast, as a 
function of dimerization 8 for impurities on the first and 
fourth dimers (a,b) and first and fifth (c), at fixed J2 = 0.5Ji. 
Also shown in (a) is the singlet-singlet separation, Ass, for 
the case of impurites on the same J2 chain, (b) shows the 
same data as (a), but magnifies the low-energy sector. 



place nearest-neighbor dimers on any link, by penalizing 
those sections where the dimer pattern preferred by the 
impurity positions lies on the weaker Ji(l — 5) bonds. 

It is immediately apparent from Fig. 7 that the evo- 
lution of the low-lying states with dimerization is not 
monotonic, and as in the ladder limit is not governed in 
magnitude or sign by a single rule. That the ladder- and 
MG-limit ground states do not evolve directly into ea- 
chother is no surprise, as a significant rearrangement of 
local dimers is necessary (Fig. 2). Similarly, the minimum 
in the singlet-triplet energy separation at 0.2 < 8 < 0.3 
arises because reducing 8 from the ladder limit increases 
the frustration between rung singlets and reduces the ef- 
fective interaction until the MG limit is reached at small 
8. 
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To investigate the nature of the intermediate dimcr- 
ization regime in Fig. 7, we have calculated the expec- 
tation values S m = (Sj-Sj+ m ) for local dimer formation 
on nearest- (m = 1) and next-nearest-neighbor (m = 2) 
sites in the chain. In a perfect singlet S\ = —3/4, and we 
find that dimer formation on the strong bonds Ji(l + <5) is 
dominant, meaning Si ~ —0.7, for all values of the dimcr- 
ization S > 0.1. The crossover from a MG pattern of flex- 
ible dimer formation occurs in the range 0.05 < S < 0.1. 
This is fully consistent with the energy spectra (Fig. 7), 
where the existence of a low-energy manifold contain- 
ing only one singlet and one triplet is established at and 
above 5 = 0.1. 




FIG. 8. Representation of dimer polarization for the 
dimerized, frustrated chain, (a) Ladder limit, inter-dimer 
coupling determined by J2 bonds: Si > between dimers, 
S2 < 0. (b) MG chain limit, inter-dimer coupling determined 
by Ji(l — 8) bonds: Si < between dimers, <S 2 > 0. 

A more detailed understanding of the effective interac- 
tions can be obtained from Si on the bonds between the 
strong dimers, and from In a ladder, S2 < because 
of the AF chain bonds, and Si > as a consequence, 
although both quantities are small [O(10~ 2 )]. In a MG 
chain with two impurities, we find that the inter-dimer 
coupling is controlled by Ji(l — 5) rather than by the 
two bonds J2, and Si < while S2 > 0. These two situa- 
tions are represented in Fig. 8. On the quantitative level, 
even for <5 ~ 0.1, the inter-dimer |<Si| and | «S*2 1 arc of or- 
der 0.1, as would be expected from the magnitude of the 
singlet-triplet energy separation (effective inter-impurity 
interaction) . 

Returning to Figs. 7(a) and (b), in the ladder limit <5 — > 
1 the ground state is a singlet (triplet) for impurities on 
the same (opposite) chains. This situation is maintained 
as 5 is lowered, until for the configuration ((l,b)(4,b)) the 
triplet crosses the singlet at <5 = 0.5 to become the ground 
state. This result has a ready interpretation in terms of 
the "free-spin" physics which governs the behavior of the 
doped ladder: when the impurities occupy the same chain 
the number of spins between them is odd, and one is a free 
spin in the ladder limit. This spin is coupled to one site of 
the neighboring dimer in the shorter inter-impurity chain 
segment by J 2 , and to the other by Ji(l — S) (Fig. 9), 
so when the latter becomes stronger than the former the 



sign of the effective interaction changes. Note that the 
other free spin does not have a Ji(l — 5) bond to the 
inter-impurity segment, and so the sign of its coupling is 
unaltered. When the impurities are on opposite chains, 
neither or both, but not only one, can have a Ji(l— S) link 
to the short segment, so that no sign-change can occur. 
The triplet ground states for dimerizations 0.25 < S < 0.5 
thus have a consistent explanation in terms of free spins 
introduced by each impurity. 

However, from the previous discussion of dimerization 
expectation values S m , the singlet regimes directly be- 
low S — 0.25 for impurity configurations ((l,b)(4,m)) 
should not be simply those of the MG limit, and this 
crossover should not mark the onset of the regime gov- 
erned by degeneracy-induced screening of the spins ex- 
plored in Fig. 5. The computed dimerization patterns 
show that for short inter- impurity segments, while dimer 
formation on the Ji(l + S) bonds remains strong, the re- 
gion around S = 0.25 does mark the crossover between 
the two types of dimer polarization. For a ((l,b)(4,m)) 
impurity configuration, where there are two dimers in the 
inter-impurity segment, dimer repolarization from a MG 
"head-to-tail" pattern [Fig. 8(b)] to a ladder-like "head- 
to-head" situation [Fig. 8(a)] accounts directly for the 
change in the spin sector of the ground state seen at 
S = 0.25 in Fig. 7(b). 




FIG. 9. Reversal of the effective interaction sign for 
odd-length inter-impurity segments at 8 = 0.5: the "free spin" 
on the right has competing couplings Ji(l — 8) and J2 to the 
neighboring singlet, which is (exponentially) weakly polarized 
by the other free spin. 

In Fig. 7(c), the behavior in the MG regime is iden- 
tical to Fig. 7(a) (see Fig. 6), and so has been cut off 
for greater clarity in the ladder limit, where the ener- 
getic separation of the singlet and triplet is now much 
smaller (Fig. 4). In the ladder regime one again sees 
singlet or triplet ground states due to the effective inter- 
action between the free spins, whose sign is selected by 
the geometry. Here it is the triplet ground state which 
becomes a singlet as (5 is reduced below 0.5, mirroring 
the behavior in Fig. 7(a,b). However, in this case both 
singlet ground states in the regime 0.25 < S < 0.5 do 
not cross directly to the singlet states of the MG limit, 
but are separated from these by a region with very weak 
triplet ground states for both "even" and "odd" impurity 
configurations, at 0.15 < <5 < 0.23. This triplet region is 
found to be robust in chains of 18 and 22 sites, which do 
not alter the short inter-impurity path, but is not present 
for the nearest accessible impurity configuration, which 
is two impurities on the first and third rungs. 
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From above, the explanation of the triplet regime is 
expected within a free-spin scenario. In contrast to the 
((l,b)(4,m)) case, the ground state for impurity configu- 
rations such as ((l,b)(5,m)) (Fig. 8), with odd numbers 
of dimers in the inter-impurity chain segment, should 
not be affected by the physics of dimcr repolarization 
since the number of inter-dimer bonds is even. Indeed 
the dimerization patterns Si and £2 are characteristic of 
the ladder polarization for dimerization values S > 0.08 
for this length of segment. In the vicinity of the triplet 
regime, the patterns are qualitatively the same on both 
sides of the crossings at 5 = 0.15 and 5 — 0.23, and 
do not indicate a difference in the nature of the ground 
state. The extremely close competition between singlet 
and triplet reflected in the very small values of their 
separation in this regime suggests that the emergence 
of a triplet ground state is a higher-order effect which 
cannot be explained using only the bond spin correla- 
tion functions S m or on-site expectation values (S z ). We 
have found a triplet ground state also for the open sys- 
tem of two "impurity" spins coupled to six spins in a 
frustrated and dimerized segment when 5 — 0.2. We 
have further verified using a 10-spin intermediate seg- 
ment that this triplet ground state is not present for 
longer inter- impurity segments with odd dimer number, 
although in this case the triplet energy approaches to 
within 2xlO~ 6 Ji of the singlet. 

To summarize this section, we have analyzed the evolu- 
tion of the free spins induced by nonmagnetic impurities 
in the spin ladder on tuning the dimerization to the MG 
limit. As frustration is increased, free-spin nature sur- 
vives down to S ~ 0.1, below which the free spins are 
screened and the low-energy sector changes rapidly to 
reflect the dominance of local dimer formation and re- 
arrangement. Although the free spins are only weakly 
bound in the ladder case, while all neighbouring spins 
are strongly bound in the MG chain, in fact the overall 
magnitude of the effective coupling is reduced towards 
zero between these limits by the presence of frustration, 
and its sign changes according to details of the impurity 
locations. 



V. PHASE DIAGRAM 

For the thermodynamic limit of infinite chains with 
random doping by a finite concentration of nonmagnetic 
impurities, the system will contain finite segments of 
all possible interjrnpurity spacings governed by a Pois- 
son distribution.Ej Our results allow us to propose the 
schematic phase diagram in Fig. 10(a). For the major- 
ity of the parameter space, the system is dominated by 
geometrical dimer formation, and for all but the low- 
est energy scales will exhibit a free spin degree of free- 
dom per impurity as these dimers are broken. The low- 
energy thermodynamic response would show the pres- 
ence of these spins and the absence of a true spin gap. 



For small dimerization and appreciable frustration, there 
is a significant regime where the impurities are screened 
and do not feature in the thermodynamic response of 
the system, which as a result may still be appropriately 
characterized as a spin liquid. For a finite region of small 
J2 and (5, we anticipate an average of half a free spin 
per impurity, simply as the result of producing effec- 
tively isolated chain segments of odd and even lengths. 
Because no true transitions are possible in the system 
at finite doping, these regions should be connected by 
crossover regimes, where we cannot exclude the possibil- 
ity of a mixed behavior, in that configuration-dependent 
local dimer formation may favor screening of some impu- 
rities while isolated spins are left around others. 
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FIG. 10. Schematic phase diagram showing the number of 
free spins per impurity for the dimerized, frustrated spin chain 
in the thermodynamic limit, (a) Finite impurity concentra- 
tion: the dashed lines do not represent true phase transitions, 
but rather denote crossovers between different regimes of be- 
havior, (b) Limit of vanishing impurity concentration. 

Fig. 10(b) shows by contrast the phase diagram which 
we conjecture for a finite number of impurities, whose 
concentration vanishes in the thermodynamic limit. We 
stress that the results S\ < and S2 > which we obtain 
for the inter-dimer correlations in the MG chain must de- 
pend on the presence of impurities, as the exactness of 
the dimerized wave function guarantees that these quan- 
tities are zero in a pure chain. In the limit of infinitely 
separated impurities, any finite dimerization 5 will then 
act to confine spin excitations, and place the system in 
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the regime of free-spin physics. This statement may be 
justified by considering the persistence of an exact dimer 
wave function, without double degeneracy, along the line 
2J2/J1 + S = 1. Thus in the somewhat hypothetical sit- 
uation of vanishing impurity concentration, the regions 
where an impurity introduces no free spins, or half of one, 
would be reduced to the S — axis [Fig. 10(b)], and the 
three regions would be separated by true transitions. 

VI. CuGe0 3 

One possible experimental realization of the dimerized, 
frustrated chain is the compound CuGe03, which is best 
known for exhibiting a spin-Peierls transition. In this 
material the dimerization is weak, and the frustration 
remains ambiguous due to competing models, both one- 
and two-dimensional, which have been used to fit the 
measured spin susceptibility. It is important to note also 
that in our model above, the dimerization parameter 6 
was taken to be rigid, meaning unaffected by the intro- 
duction of a nonmagnetic impurity, which corresponds 
to the limit of very high phonon frequencies. A further 
point of contention in CuGeOs is the rigidity of the lat- 
tice, or the extent to which the local dimexization pattern 
is in fact relaxed around impurity sites £3 

The above uncertainties notwithstanding, the most di- 
rect experimental measure of the number of free spins in 
a system is the Curie- Weiss tail in the low-temperature 
susceptibility. This has been reported to scale with the 
impurity concentration in some ladder systems and in 
the spin-Peierls phase of CuGe03, implying the creation 
of one free spin per impurity. However, highly doped 
samples of CuGeO-3 , which do not undergo a spin-Peierls 
transition, do not show a significant Curie- Weiss tail, and 
in any case nothing .which may be scaled with the im- 
purity concentrationO We suggest that the presence of 
frustration, which is in addition likely to be stronger in 
the undimerized phase than in the spin-Peierls phase, is 
a viable candidate to explain this situation. 

Our analysis implies a lower bound on this frustration 
of J2/J1 > 0.25, a conclusion which cannot be drawn 
from the susceptibility alone. A further qualitative de- 
duction is that frustrating interactions are much more 
significant than interchain interactions, which also act 
to stabilize the dimer pattern and free-spin physics, in 
this doping regime. Finally, we note that an impurity- 
screening mechanism due only to lattice relaxation would 
result on average in screening of only half of the impuri- 
ties. 



VII. CONCLUSION 

In summary, we have analyzed the behavior on doping 
by nonmagnetic impurities of the dimerized, frustrated 
Heisenberg spin chain. We find that when frustration is 



predominant it can act to screen the impurity-induced 
spins, even in systems with a spin gap. This is in con- 
trast to the cases of the spin ladder and dimerized chain, 
where one impurity acts to introduce one free spin. In 
the general case of a chain with dimerization and frus- 
tration, free-spin physics around the impurities breaks 
down when the dimerization is reduced below 8 <~ 0.1, to 
be replaced by a screening of the impurities which arises 
from the near-degeneracy of the J\ bonds. The sign and 
magnitude of the effective interactions show a complex 
but systematic evolution between the limiting cases. For 
long chains with finite impurity concentrations we deduce 
a phase diagram with crossovers, but no true phase tran- 
sitions, separating regions where each dopant introduces 
0, 1/2, or 1 free-spin degree of freedom. 

More generally, we may conclude that nonmagnetic 
impurities in quantum magnets give rise to a variety of 
subtle phenomena. The essential feature of the undoped 
MG chain which leads to screening of the impurity sites 
is the presence of two degenerate singlet ground states 
with a finite gap to the first triplet. We may expect sim- 
ilar screening effects, leading to an absence of free-spin 
degrees of freedom, to operate around nonmagnetic im- 
purities in other frustrated spin systems with analogous 
properties. A notable example in this category is the S 
= 1/2 Heisenberg model on the Kagome lattice, where 
the ground-state degeneracy takes the form of a contin- 
uum of low-lying singlets in the singlet-triplet gap.Ell Be- 
cause a resonant valence-bond (RVB) formulation based 
on singlet dimer coverings has been shown to give a good 
description of these low-lying singlets,Ed the same general 
screening mechanism, which is based on rearrangement 
of singlets around the impurities, is expected to apply 
also to this case. 
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